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Abstract 



The smoothed correlation function for the eigenvalues of large hermitian 
matrices, derived recently by Brezin and Zee [Nucl. Phys. B402 (1993) 613], 
is generalized to all random-matrix ensembles of Wigner-Dyson type. 



O ■ I- INTRODUCTION 

o ■ 

^ ■ A basic problem in random-matrix theory is to compute the correlation of the eigenvalue 

^ . density at two points in the spectrum, from the Wigner-Dyson probability distribution of 
\ the eigenvalues [||]. The correlation is a manifestation of the level repulsion resulting from 
the jacobian Yli^j \\i — Aj|^, which is associated with the transformation from the space 
of N X N hermitian matrices to the smaller space of N eigenvalues Ai, A2, . . . Aat. [The 
■ power (3 depends on whether the matrix elements are real {(3 = 1, orthogonal ensemble), 
O ■ complex [jS = 2, unitary ensemble), or quaternion real {jS = 4, symplectic ensemble).] The 



Wigner-Dyson probability distribution 



X- P{{Xn}) = Z-' exp[-l3W{{X^})], (la) 

H ; 

with Z a normalization constant and 

W{{Xn}) = - E In |A^ - A, I + 5: V{X,), (lb) 

i<j i 

describes an ensemble where all eigenvalue correlations are due to the jacobian. The poten- 
tial V{X) determines the mean density p(A) of the eigenvalues, which is non-zero in some 
interval (a, b). 

In many applications of random-matrix theory, it is sufficient to know the eigenvalue 
correlations in the bulk of the spectrum, far from the end points at a and b. In some 
applications, however, the presence of an edge in the spectrum is an essential part of the 
problem, and its effect on the spectral correlations can not be ignored. The application to 
"universal conductance fluctuations" in mesoscopic conductors is one example The 
application to random surfaces and two-dimensional quantum gravity is another example 
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Recently, Brezin and Zee |P,|TO[ reported a remarkably simple result for the two-level 
cluster function 



T,{\ ^) = - ( ^ 5(A - A,)5(/x - A,) ) + p(A)p(/i), (2) 

which included the effects of an upper and lower bound on the spectrum. (Here (■ ■ ■) denotes 
an average with distribution ([l|), and p(A) = (X^i ^(A — Aj)) is the mean eigenvalue density.) 
For ^ 1, the correlation function oscillates rapidly on the scale of the spectral band 
width {a,b). These oscillations are irrelevant when integrating over the spectrum, so that 
in the large- limit it is sufficient to know the smoothed correlation function. Brezin and 
Zee considered the unitary ensemble {(3 = 2), with V{X) = I]fc=i CfeA^*^ an even polynomial 
function of A, so that a = —b. (The case a ^ h can then be obtained by translation of 
the entire spectrum.) Using the method of orthogonal polynomials they computed the 
smoothed correlation function, with the result 

\^ _}__}: - ^/^ rr>\ 

' ^' 27r2 (A - p)2 [(a2 - \^){a^ - p2)]i/2 ' 

The purpose of the present paper is to show how Eq. (|^) can be generalized to arbitrary 
(non-polynomial, non-even) potentials ^(A), and to all three symmetry classes {(3 = 1,2,4). 
This universality is achieved by a functional derivative method , which provides a powerful 
alternative to the classical method of orthogonal polynomials. In Ref. |^ we applied this 
method to the case a = 0, 6 — oo of a single spectral boundary. Here we extend the analysis 
to include a finite upper and lower bound on the spectrum. 



II. METHOD OF FUNCTIONAL DERIVATIVES 

We consider the two-point correlation function 

K,{X, /i) = - ^(^ - ^^n^^ - ^j)) + PWpif^) (4) 

(note the unrestricted sum over i and j), which is related to the two-level cluster function 
(0) by 

K^iX, fi) = T2(A, p) - p(A)5(A - p). (5) 

For A 7^ p the two correlation functions coincide, so that we can compare with Eq. (^. 
We prefer to work with K2 instead of T2 for a technical reason: Smoothing, in combination 
with the large- limit, introduces a spurious non-integrable singularity in T2(A, p) at A = p, 
while K2{X,fi) remains integrable.0 



^The distinction between T2(A, /u) and i^2(A,p) was not made explicitly in Ref. Q, because only 
the case A / was considered. 
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Our analysis is based on the exact relation between the two-point correlation function 
K2{X, fi) and the functional derivative of the eigenvalue density p(A) with respect to the 
potential V{fi), 

The smoothed correlator is obtained by evaluating the functional derivative using the asymp- 
totic (A^ oo) integral relation between V and p, 

r fdp{^ = ^ViX), a<X<b. (7) 
Ja A — /i aX 

(The symbol V denotes the principal value of the integral.) Corrections to Eq. are 
smaller by an order for /3 = 1 or 4, and by an order A^~^ for j3 = 2 [O]. Variation of 
Eq. gives 

s,m. ,,i^^v!V-^=^jvw, (8) 

X — b X — a Ja X — n dX 



with the constraint 



dA 5p{X) = (9) 



(since the variation of p is to be carried out at constant A^). The end point a is either a 
fixed boundary, in which case 6a = 0, or a free boundary, in which case p{a) = 0. Similarly, 
either 56 = or p{b) = 0. We conclude that we may disregard the first two terms in Eq. (^, 
containing the variation of the end points. What remains is the singular integral equation 

-p ['dp = -^SViX), a<X<b, (10) 

J a X — p dX 

which we need to invert in order to obtain the functional derivative 6p/6V. 
The general solution to Eq. ([To|) is [|l^ 

^'^^^ = ^ [(A-a)(5-A)]V^ [^-^Ja A^^;^ ^^^^^^ j " ^''^ 

The coefficient C is determined by 

C = n fdX5p{X). (12) 



In view of Eq. (H), we have C = 0. Combination of Eqs. (H) and ([TT|) yields the two-point 
correlation function 
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(We have substituted the representation Vx~^ = (d/dx) In |x| for the principal value.) 

The two-point correlation function (0) has an integrable singularity for A = /i. For 
A 7^ one can carry out the differentiations, with the result 



^2(A,/i) 



1 1 l{a + b){X + fi) - ab - Xfi 

^ (A - /i)2 [(A - a){b - A)(^ - a){b - /i)]i/2 



(14) 



For a = —b and /5 = 2 we recover the formula (|^) of Ref. for even polynomial potentials 
in the unitary ensemble. For a = and 6 ^ oo we recover the correlator of Ref. 0, 



^2(A,/i) 



d d 



jSn"^ dX djj, 
1 



In 



v/A - v^Ai 



^A + 

(A - ^)-2(A + /i)(A^)-^/2 if A ^ ^, 



(15) 



for the case of a single spectral edge. 

An important application of the smoothed two-point correlation function is to compute 
the large- limit of the variance VarA = {A^) — {A)"^ of a linear statistic A = J2n=i '^i'^n) 
on the eigenvalues, by means of the relationship 



pb rb 

VarA = — dX dfia{X)a{fi)K2{X, fi). 

J a J a 



Substituting Eq. (|13|) we obtain, upon partial integration. 



VarA 



V dX dfi 



{jj, — a){b — jj) 
(A-a)(6- A) 



1/2 



a(A) d 



X — fi dfi 



(16) 



(17) 



Note that here it is essential to work with the expression (^) for K2{X,fi), which is inte- 
grable, and that one can not use the expression (|T^), which has a spurious non- integrable 
singularity at A = /x. The formula (p!7[) is the generalization to a spectrum bounded from 
above and below of previous formulas by Dyson and Mehta [|T^ (for an unbounded spectrum) 
and by the author (for a spectrum bounded from below). 



III. CONCLUSION 

The result (|13]) for the smoothed two-point correlation function in the large- iV limit holds 
for all random-matrix ensembles of Wigner-Dyson type, i.e. with a probability distribution 
of the general form (|l]). The form of the eigenvalue potential V{X) is irrelevant. It is also 
irrelevant whether the end point at A = a (or at b) is a fixed or a free boundary. This is 
remarkable, because the eigenvalue density behaves entirely different in the two cases: At 
a fixed boundary p(A) diverges as (A — a)"^''^, while at a free boundary p(A) vanishes as 
(A — a)^/2. Both cases are of interest for applications: The spectrum considered in Ref. 0, 
in connection with two-dimensional gravity, has free boundaries; The spectrum considered 
in Ref. ||^, in connection with mesoscopic conductors, has a fixed boundary. 

While the form of the eigenvalue potential is irrelevant, the form of the eigenvalue inter- 
action does matter. Consider an eigenvalue distribution function of the form (|ip, but with a 
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non-logarithmic eigenvalue interaction u{X, jj,) 7^ In |A — Such a distribution describes the 
energy level statistics of disordered metal particles [14|, and the statistics of transmission 
eigenvalues in disordered metal wires |]T5|. The analysis of Sec. 2 carries over to this case, 
but the integral kernel (A — in Eqs. (|^) and (|l3) has to be replaced by the kernel du/dX. 
The two-point correlation function now equals 1/(3 times the inverse of this integral kernel, 
and differs from the result (0) for a logarithmic interaction. 

So far we have only considered the two-point correlation function K2 = (3~^5p/5V and 
the closely related two- level cluster function T2. Brezin and Zee also computed the three- 
and four-level cluster functions, and found that they vanish identically upon smoothing. 
The linearity of the relation (|^) between p and V implies in fact that this holds for all 
higher-order cluster functions. This argument is equivalent to Politzer's proof |jl6| that any 
linear statistic on the eigenvalues has a gaussian distribution in the large- iV limit. 
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